Anomalous optical and electronic properties of dense sodium 
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Abstract 

Based on ab initio density-functional-theory using generalized gradient approximation, we systematically study the optical 
and electronic properties of the insulating dense sodium phase (Na-hp4) reported recently [Ma et al., Nature 458, 182 (2009)]. 
The structure is found optically anisotropic and transparent to visible light, which can be well interpreted using its electronic 
band structure and angular moment decomposed density of states. Through the bader analysis of Na-hp4 at different pressures, 
we conclude that ionicity exists in the structure and becomes stronger with increasing pressure. In addition, the absorption 
spectra in the energy range from 1.4 to 2.4 eV are compared with recent experimental results and found good agreement. It is 
found that the deep-lying valence electrons participate in the interband transition. 

PACS numbers: 71.20.-b, 71.20.Dg, 78.20.Ci, 78.20.-e 
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I. INTRODUCTION 



The electronic structure and property of alkali met- 
als under pressure have drawn extensive attentions re- 
cently. At low pressures they all possess simple body- 
centered-cubic (bcc) or face-centered-cubic (fee) struc- 
ture, for which the tightly bound core electrons remain 
largely unaltered with respect to those in the free atoms 1 . 
However, the application of pressure may cause the core 
states to come into play and thus results in unexpected 
varieties of complex phenomena. Sodium (Na) is unique 
among alkali metals due to the occupied electronic states. 
Under pressure, Na is apt to adopting low-symmetry 
structures^ 3 - and even incommensurate phased. The high 
pressure transition sequence of Na has been recognized as 
follows: Na in bcc structure transforms to fee structure at 
65 GPa 6 and then to a more complex cI16 structure (16- 
atom bcc) at 105 GPa 5 -! 7 -; The orthorhombic oP8 phase 
is observed between 118 to 125 GPa 3 ; At 125 GPa, oP8 
transforms to an incommensurate phase tI19^. Though 
the darkening of Na at the transition to til 9 was ob- 
served, the phase is still notably metalio^. Until recently, 
Ma et al. reported a new phase of Na transparent to 
visible light at pressure about 200 GPa s , which suggests 
that the metal has become an insulator or at least semi- 
conductor. First-principles evolutionary methodology for 
crystal structure prediction identified the new phase as 
wide bandgap dielectric with a six-coordinated, highly 
distorted double-hexagonal close-packed structure, de- 
noting as hp4 8 . This finding may stimulate the study 
of transition mechanisms from metal to nonmetal un- 
der high pressures and promote the development of high 
pressure theory as well. It is encouraged to study the 
electronic properties of the insulating phase Na-hp4. 



Optical measurement is a powerful tool for studying 
the electronic structures of semiconductor and metal. 
The optical absorption and reflective spectra of Na 
have been studied extensively both experimentally and 
theoretically 4 ^— . In particular, Lazicki et al. have 
recently studied the optical properties of Na in bcc, 
fee, cI16, oP8 and tI19 phases comparatively using 
synchrotron infrared spectroscopy and first-principles 
method within DFT 1 - 2 . It was found that Na gradu- 
ally loses its Drude-like luster under compression, and 
tends towards less metallic behavior. As for the insulat- 
ing phase Na-hp4, though spectroscopic data have been 
obtained, the main optical spectra are still scarce except 
experimentally measured absorption spectra in a small 
energy range. It is indispensable to study the optical 
properties of Na-hp4 systematically both for interpreta- 
tion of the mechanism of metal transition to insulator 
and for further insight into the electronic structure of 
Na-hp4. As dielectric material, exploration of dielectric 
constants of Na-hp4 may be useful from the viewpoint 
of technological applications. Besides, it is also expected 
that the calculated results can be as useful reference for 
future experimental studies. 

In the present study, the electronic structure and opti- 
cal properties of Na-hp4 are calculated based on ab initio 
DFT using generalized gradient approximation (GGA). 
These calculations demonstrate that Na-hp4 is transpar- 
ent to visible light anisotropically. The p-d hybridization 
of valence electrons and strong interstitial charge concen- 
trations are the main character of Na-hp4. Through DOS 
and Bader analysis of Na-hp4, we conclude that electrons 
transfer between the two inequivalent atoms in the unit 
cell and ionicity exists between them. And also the ion- 
icity becomes stronger with increasing pressure. The rest 
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of this paper is organized as follows. In the next section, 
the computational details regarding the methods used in 
our calculations of electronic structure and optical prop- 
erties are described. In Sec. Ill, the calculated results 
are discussed and compared with experimental data£. Fi- 
nally, we close our paper with a summary of our main 
results. 



II. COMPUTATIONAL DETAILS 

Our electronic and optical calculations are performed 
using a plane-wave implementation of DFT within the 
GGA of Perdew-Burke-Ernzerhof (PBE) formalism^, as 
implemented in the VASP code (Vienna ab initio simu- 
lation program)*™. The all-electron projector augmented 
wave (PAW) metho d 15 i 16 are adopted, with the PAW po- 
tential treating 2s 2 2p 6 3s 1 as valence electrons. It is im- 
plemented by incorporating Is 2 state into an effective 
frozen core. The PAW potential provides a significant 
advantage over the ultrasoft pseudopotential in applying 
to the optical property calculation since the all electron 
PAW treatment avoids the complication of calculating an 
awkward correction term related to the nonlocal pseu- 
dopotential operato r 17 i 18 . We use the plane-wave kinetic 
energy cutoff of 1250 eV, which is shown to give excellent 
convergence of the total energies. For structure relax- 
ation and electronic properties calculations, 12x12x10 
Monkhorst-Pack grid are used, while for optical property 
study, the T-centered 18x18x12 k grids and 72 bands are 
adopted yielding dynamical dielectric constants within 
0.002 (compare Table II). 

The response of a system to an external electromag- 
netic field with a small wave vector can be described with 
the complex dielectric function. We consider only elec- 
tronic excitations because only the dielectric function for 
frequencies well above those of the phonons is of interest 
to us. For these, two different methods are used to deter- 
mine the dielectric constants using different approxima- 
tions. On one hand, the linear response theory (density 
functional perturbation theory) is applied; on the other 
hand, a summation over conduction band states is used. 
For the former, only the static ion-clamped dielectric ma- 
trix can be obtained without summation of conduction 
bands required. While the latter can be used to calcu- 
late the frequency-dependent dynamical dielectric con- 
stants when Kohn-Sham eigenvalues and eigenfunctions 
are obtained. Moreover, due to the double- hexagonal 
close-packed structure of sodium, the dielectric function 
is an anisotropic tensor. By an appropriate choice of the 
principal axes, we can diagonalize it and restrict our con- 
siderations to the diagonal matrix elements. The inter- 
band contribution to the imaginary part of the dielectric 
function is calculated by summing transitions from oc- 
cupied to unoccupied states (with fixed k vectors) over 
the Brillouin zone, weighted with the appropriate matrix 
element giving the probability for the transition. To be 
specific, the longitudinal expression for the components 



e^g (w) are given by 
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where the indices c and v refer to conduction and valence 
band states respectively, and w c k is the cell periodic part 
of the wavef unctions. The vector e Q are the unit vectors 
for the three Cartesian directions. SI is the volume of the 
primitive cell. While the k— point weights are defined 
such that they sum to 1, the Fermi weights equal to 1 for 
occupied and zero for unoccupied states. The factor 2 
before the weights accounts for the fact that we consider 
a spin-degenerate system. The real part of the dielectric 
tensor e (w) is obtained by the usual Kramers-Kronig 
transformation: 
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where P stands for the principal value of the integral. 
The knowledge of both the real and imaginary parts of 
dielectric tensor allows the calculation of important op- 
tical constants. In this paper we present and analyze 
the reflectivity R (uj), the absorption coefficient I (ui), the 
electron energy-loss spectrum L (w), as well as the refrac- 
tive index n and extinction coefficient k. The reflectivity 
spectra are derived from the Fresnel's formula: 
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The absorption coefficient / (w) and the electron energy- 
loss spectrum L (w) are expressed explicitly as, 
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The refractive index n and the extinction coefficient k 
are obtained by the following expressions: 
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Figure 1: (Color online) (a). The unit cell of Na-hp4 struc- 
ture, x, y, and z are the crystallographic axes used for prop- 
erty studies. Nal and Na2 are the two inequivalent atoms in 
the unit cell, (b) Top view of the Na-hp4 structure. A, B and 
C correspond to the different layers shown in (a). 




The optical conductivity then follows immediately from 
the imaginary part of the dielectric constants, 



Figure 2: (Color online) Band structure of Na-hp4 at 320 
GPa. The Fermi energy was set at the top of the valence 
band. 



^( W ) = - £ $( W ). (8) 



III. RESULTS AND DISCUSSION 

A. Atomic and electronic structures of Na-hp4 

The crystal structure of sodium in hp4 phase have 
been studied experimentally using x-ray diffraction and 
reflection 8 . The insulating phase considered here is a 
highly distorted double-hexagonal close-packed structure 
(D_2h space symmetry). It belongs to the space group 
of P&s/mmc and contains four atoms per cell. There 
are two inequivalent atomic positions in unit cell, Nal 
located at the 2a site (0.0,0.0,0.0) and Na2 at 2d site 
( 2 /3, !/3, All the crystal parameters are optimized. 
Atoms in this dense structure are six-coordinated. Note 
that the stacking of close-packed layers of Na atoms is 
CACBCACB. . . as in any d.h.c.p structure. Top view 
along the c-axis within 2a x 2b plane of Na-hp4 is show 
in the Fig. 1(b), in which A, B, and C denote the atoms 
occupying the A, B, and C layers in the Fig. 1(a). 

Electronic structure is first calculated since the op- 
tical spectra are calculated from interband transition. 
It should be noted that theoretical predicted transition 
pressure of Na-hp4 is about 260 GPa, which is higher 
than experimental reported value 200 GPa 8 . The dis- 
crepancy may result from the well-known overstabiliza- 
tion of metallic states by DFT calculations. For this rea- 
son, we calculate the electronic and optical properties of 
Na-hp4 at 320 GPa, at which Na-hp4 exists stably. The 
calculated band structure of hp4 phase at 320 GPa in 
the high-symmetry directions in the Brillouin zone are 
shown in Fig. 2. The energy scale is in eV, and the 
valence band maximum is set to be the Fermi energy, in- 
dicating that the bands below Ef are all occupied. From 



this figure it is clear that a direct band gap appears be- 
tween the top-most valence band and the bottom-most 
conduction band at T. Our calculated value of the direct 
band gap is about 1.75 eV, which should be smaller than 
the true value due to the discontinuity of the exchange- 
correlation potential, while consistent with other theo- 
retical results 8 . The bands have large dispersion. There 
are only two valence bands in the valence region shown 
in Fig. 2, the bottom one mainly arising from 3s orbitals 
and the above one from the hybridization between 3p 
and 3d orbitals. These two valence bands are degenerate 
along A — L — H — A directions. 

Our calculated total density of states (DOS) along with 
site-projected and angular momentum decomposed DOS 
of hp4 phase are shown in Fig. 3. The valence maximum 
is set to be zero on the energy scale. The most intriguing 
feature of the DOS is the emergence of 3d state. Un- 
der the compression of 320 GPa, the overall density is 
increased to about six times as large as the atmospheric 
density of Na. The nearest Na-Na distance is 1.878 A, 
while the 3s and 2p orbital radii are 1.71 and 0.28 A, 
ionic radius of Na + is 1.02 A. That means strong valence- 
core and core-core overlap exist between the neighboring 
atoms. Due to the exclusionary effects, the Bloch states 
with the angular momentum components of 3s and 3p 
are largely excluded from the core region by the core 
states with the same angular momentum components. 
However, the more asymmetric states which have lower 
kinetic energy is apt to sampling the full nuclear poten- 
tial. Thus 3d states drops below the Fermi energy and 
hybridize with ip states, which induce strong electron 
localization as seen in valence charge density shown in 
Fig. 4. Figure 3(b) indicates that the DOS for p x +p y is 
different from that for p z , which predicts the anisotropy 
of Na-hp4. As shown in Fig. 3(c), the DOS peaks for 
Nal and Na2 are largely different, which implies evident 
charge transfer from Nal to Na2 and thus weak ionicity 



3 



— total s p 
" _p _d / 




i 

(a) A 

A 




(b) ' 


— i — i — i — i — i — i — 

. _Na1_total 

— Na2_total L 







-8 -6 -4 -2 2 4 



O 0.16 



'■_ — Na1_s 
; — Na2_s 




: — Na1 p ft 
: — Na2 p / 




■ — Na1_d 

■ — Na2 d 





i i i Z^J ■ r i , i , 

-8 -4 4 8 12 

Energy (eV) 



Figure 3: (Color online) Calculated total density of states, 
site-projected and angular momentum decomposed density of 
states of Na-hp4 at 320 GPa. The zero of the energy scale 
shows the position of the Fermi level. 



exists between them. In the valence band, the s state 
covers larger energy range than the p and d states. From 
Fig. 3(d)-3(f), the s, p and d states decrease to zero at 
7.4, 5.8 and 3.1 eV below Ef, while begin to emerge at 
2.3, 1.8 and 2.2 eV above Ef, respectively. The contribu- 
tions from p and d states imply pd and pd 2 hybridizations 
for Nal and Na2, respectively. The 3d contribution rises 
basically with increasing binding energy in conduction 
band, while the s state contribution begins to decrease 
with energy at about 3 eV above Ef. 

In order to clarify the cause of the insulating phase, 
we present in Fig. 4 the valence charge density in the 
(110) plane of Na-hp4. Density isocontours are drawn 
at intervals of 0.05 eA color-coded from (blue) to 
1 eA (red). It shows that the charge concentration 
is very strong in interstitial regions, while minimal near 
and between the ions, which is resulted by the exclusion- 
ary effect mentioned above. The transition to the insu- 
lating phase (Na-hp4) is eventually culminated. From 
this figure, charge transfer from Nal to Na2 can also be 
concluded consistently. To describe the ionic character 
quantitatively and more clearly, we perform the Bader 
analysis of Na-hp4 at different pressures. Bader analysis 
is a well established analysis tool for studying the topol- 
ogy of the electron density and suitable for discussing 
the ionicity of a materia l 19 ! 20 . The charge Qb enclosed 
within the Bader volume (Vb) is a good approximation to 
the total electronic charge of an atom. In present study 
128 x 128 x 192 charge density grids are used, which give 
excellent precision of the effective charge (< 0.01%). The 
calculated results are listed in Table I. From Table I, the 
ionic character is obviously shown through a charge flux 




Figure 4: (Color online) Valence charge density for Na-hp4 in 
(110) plane at 320 GPa. Contour lines are drawn from 0.0 to 
1.0 with interval of 0.05 eA" 3 . 

Table I: Calculated charge and volumes according to Bader 
partitioning as well as Nal-Na2 distances at different pres- 
sures. 



Pressure Qs (Nal) Q s (Na2) F s (Na2) Vs(Na2) Nal-Na2 



(GPa) 
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(A-) 


(A) 


260 


8.4024 


9.5976 


5.0458 


9.1543 


0.9222 


280 


8.3950 


9.6050 


4.8979 


8.8720 


0.9191 
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8.3875 


9.6125 


4.4607 


8.6191 


0.9067 


320 


8.3828 


9.6172 


4.6468 


8.3879 


0.8991 


340 


8.3783 


9.6217 


4.5404 


8.1746 


0.8920 



(about 0.6 electrons per Na atom) from Nal to Na2. Be- 
sides, the valence of charge transfer increases with in- 
creasing pressure, which implies that the ionicity become 
stronger under compression. The main reason lies in that 
the energy level for Na2 is lower than that for Nal. With 
increasing pressure, the energy level difference becomes 
larger. 



B. Optical properties of Na-hp4 

The interband optical functions of Na-hp4 at 320 GPa 
calculated using expressions (1) and (2) are shown in Fig. 
5(a) and 5(b), in which x, y and z correspond to the 
three principle directions E||x, E||y, and E||z, regard- 
ing to the electric field vector polarized along the crys- 
tallographic axes x, y and z respectively. It is evident 
that Na-hp4 is optically anisotropic: E||x and E||y spec- 
tra are the same, while E||z spectrum is different. The 
macroscopic dielectric constants £oo are calculated using 
different methods and approximations, as shown in Table 
II. A T-centered grid with 18xl8xl2/c points yields 
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Table II: Ion-clamped static macroscopic dielectric constants 
£00 calculated using density functional perturbation theory 
and PAW method for various fc— point sets: V indicates a 
grid centered at F, whereas Monkhorst-Pack (MP) grids do 
not include the T point, values for ej^c neglect the local 
field effects, and are obtained using linear response theory. 
£ cond are vames obtained by summation over conduction band 
states. 

Na-hp4 N k (IBZ) e^ c e cond e^r £ cond 
E||x (y) E||z 
(18xl8xl2)r 637 5.785 5.457 11.116 10.761 
(21 x 21 x 14) T 968 5.785 5.458 11.055 10.701 
(24 x 24 xl6)T 1413 5.786 5.458 11.028 10.674 
(18 x 18 x 12) MP 1512 5.786 5.459 10.900 10.546 
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Figure 5: (Color online) (a) Real part and (b) imaginary part 
of the dielectric function, and (c) absorption spectrum I (to) 
along the principle axes at 320 GPa. (d) Absorption spectrum 
I (oi) in the energy range from 1.4 eV to 2.4 eV at 200 GPa. 
In (d) the red line is theoretical value, while the black line 
represents the experimental data reported in Ref. [8]. 



dielectric constants with good convergence within 0.001 
for E||x (y) and 0.09 for E||z. The value of e^ c can 
be compared with experimental data. Obviously, Eqo is 
5.785 and 11.028 in E||x (y) and E||z directions, respec- 
tively. The obtained values can be as important reference 
for further experimental measurements of Na-hp4. 

Our calculated optical spectra cover the photon energy 
range of 30 eV. In order to show the details clearly, we 
use logarithm longitudinal scale for the imaginary part 
of dielectric constant Im (e) and absorption coefficient 
/ (u). Im (e) shows mainly five peaks both in E||z direc- 
tion and in E||x (y) directions. In E||z direction, these 
peaks correspond to the photon energy of 2.54, 4.31, 
13.37, 17.7 and 24.64 eV, while in E||x (y) directions, 
4.31, 10.54, 16.34, 20.8 and 28.8 eV, respectively. The 
peak at 2.54 eV originates from the transition from 3p 
to 3c? corresponding to the topmost valence band to bot- 
tommost conduction band, while the peak at 4.31 eV is 
attributed to the transition from 3s to 3p. At 4.31 eV 
and 13.37 eV in E||z direction and 4.31 eV and 16.34 eV 
in E||x (y) directions, the sign of Re (e) changes, which 
predict the anisotropy of the polarized Re (e) spectra. In 
the range around 13.37 eV and 16.34 eV, Im (e)is very 
low and thus the conditions for a plasma resonance are 
satisfied^. Hence, there are strong resonance maxima 
at these energies in the calculated energy-loss spectra 
L (uj) for these polarizations, as can be seen from Fig. 
6(c). The other seven high-energy peaks can all be at- 
tributed to transitions from the top two valence bands to 
the upper conduction bands, without deep-lying valence 
orbitals participating in the interband transitions. These 
also explain the origin of the peak structure in the ab- 
sorption coefficient / (u) and reflectivity R (w) spectra. 
The absorption spectra remains nearly zero in the visible 
photon energy, which implies that Na-hp4 is transparent 
to visible light. As seen from Fig. 5(c), Na-hp4 remains 
transparent to light with higher frequency in E||x (y) di- 



rection than in E||z direction. In Fig. 5(d), we present 
I (10) in the photon energy range from 1.4 eV to 2.4 eV at 
200 GPa along with the previously reported experimental 
measurement^. Here we assume the sample thickness is 
8 /im. As all know, the band gap can be deduced by ex- 
trapolating the absorption to zero. It is concluded that 
the results from theoretical calculation and experiment 
agree well except that the theoretically predicted band 
gap is smaller than the experiment, the origin of which 
is the local-density approximation. In addition, it should 
be noted that the present calculation only pertains to 
the electronic response, and does not include the effects 
of lattice vibrations that dominate the experimental ab- 
sorption spectrum in the low frequency region. Thus we 
can not derive the complete agreement with experimental 
results. 

In Fig. 6, we show the main optical spectra includ- 
ing optical conductivity, reflective spectrum, energy-loss 
spectrum, extinction and refractive indices in a photon 
energy width of 30 eV. The peak of optical conductivity 
in z direction is almost twice as large as that in x and y 
directions. The reflectivity spectrum is distinct from the 
drude-type spectrum, which means Na-hp4 is no longer 
metallic. The reflectivity is small at visible optical en- 
ergy, which ensures the transparency of Na-hp4. During 
the range about 5 eV to 10 eV, the reflective index re- 
mains almost constant as 60% in E||z direction and 20% 
in E||x (y) direction. Refractive indices of a crystal are 
closely related to electronic polarizability of ions and the 
local field inside the crystal. Another interesting point 
connected with the refractive indices is the electro-optic 
effect. We present the refractive indices n (u) in Fig. 
6(d). The low- frequency refractive index along E||x (y) 
and E||z direction are calculated as n x ^ = 2.38 and 
n z = 3.26 respectively. It is clear that the light polarized 
parallel to the z axis is more refracted than that with 
polarization along the x and y axes. This indicates the 
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Figure 6: (Color online) (a) Calculated optical conductivity 
spectrum a (u), (b) reflectivity spectrum R(oj), (c) energy- 
loss spectrum L(uj), (d) refractive index n(uu), and (e) ex- 
tinction coefficient k (u) along the principal axes at 320 GPa. 



ent insulating phase Na-hp4 using the first-principles 
method. The insulating mechanism that exclusionary ef- 
fect induces the valence electrons being excluded from the 
core region is clearly demonstrated from the calculated 




12 16 
Energy (eV) 



Figure 7: (Color online) Calculated effective number of elec- 
trons rieff participating in the interband optical transitions 
along the principle axes at 320 GPa. 



large optical anisotropy in Na-hp4. 

Since the optical properties of Na-hp4 are anisotropic, 
it is particularly interesting to calculate the effective 
number of valence electrons per Na atom n e ff contribut- 
ing to the optical properties in each direction. This can 
be accomplished by means of the sum rule^: 



2m 
Ne 2 h 2 



E (Im (e)) dE, 



(9) 



where E m denotes the upper limit of integration, m is 
the electron mass, and N stands for the electron density. 
The effective number of valence electrons per Na atom 
participating in the interband transitions along the crys- 
tallographic directions at 320 GPa is shown in Fig. 7 as 
a function of E m . One can see that n e ff remains lager be- 
fore 13 eV and smaller after 13 eV in E||z direction than 
in E||x (y) direction. This is connected with the larger 
reflectivity feature observed in E||z direction at low en- 
ergy range than in E||x (y) direction. In addition, the 
n e ff doesn't reach a saturation value up to 30 eV. This 
shows that deep-lying valence orbitals do participate in 
the interband transition. 



band structure, angular moment decomposed density of 
states and valence charge density. The dielectric func- 
tions and all main optical spectra indicate that Na-hp4 
is optically anisotropic. The macroscopic dielectric con- 
stants have been obtained and are expected to be 
useful for future reference. Especially, our calculated ab- 
sorption spectrum in the energy range from 1.4 to 2.4 
eV agrees well with the experimental observation. Fi- 
nally, we have shown that during the interband transi- 
tion, deep-lying valence electrons of Na-hp4 also partici- 
pate. 
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